Pseudospin magnetism in graphene 
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We predict that neutral graphene bilayers are pseudospin magnets in which the charge density- 
contribution from each valley and spin spontaneously shifts to one of the two layers. The band 
structure of this system is characterized by a momentum-space vortex which is responsible for 
unusual competition between band and kinetic energies leading to symmetry breaking in the vortex 
core. We discuss the possibility of realizing a pseudospin version of ferromagnetic metal spintronics 
in graphene bilayers based on hysteresis associated with this broken symmetry. 



Introduction — The ground state of an interacting elec- 
tron system flows from subtle compromises between band 
and interaction energy minimization. Because of the 
Pauli blocking effects which underlie Fermi liquid the- 
ory however, the consequences of interactions are nor- 
mally only quantitative [l[ unless symmetries are bro- 
ken. Recent progress Q in the isolation of single and 
double sheets of carbon atoms (graphene sheets) has pre- 
sented researchers with a new type of interacting elec- 
tron system whose properties are now being actively ex- 
plored [H, 0, HI, both theoretically and experimentally. 
In this Letter we argue that band energy minimization is 
exceptionally frustrating to interactions in graphene bi- 
layers, and predict that broken symmetry states in which 
charge shifts spontaneously from one layer to the other 
occur as a consequence. 

Graphene bilayers with Bernal stacking have one 
low-energy site per unit cell in each layer. When 
the layer degree of freedom is described as a pseu- 
dospin, the continuum limit of the 7r-orbital band 
Hamiltonian corresponds @, 0] to a pseudospin field 
Bba„d = [h 2 k 2 /(2m*)} (cos(20 fe ),sin(20 fc ),O), where 
4>k = arctaxi (k y /k x ), m* — 7i/(2v|), 71 is the interlayer 
tunneling amplitude, and vf is the electron velocity at 
the Fermi energy in an isolated neutral graphene sheet. 
When interactions are neglected the ground state of a 
neutral bilayer has a full valence band of pseudospinors 
aligned at each k with this pseudospin field, forming the 
momentum-space vortex. The vortex exacts a large in- 
teraction energy penalty because of its rapid pseudospin- 
orientation variation. We propose that, like its real-space 
counterpart [8], the momentum-space vortex sidesteps 
this energy cost by forming a vortex core in which the 
pseudospin orientation is out of plane in either the z or 
— z direction, as illustrated in Fig. [TJ The momentum- 
space vortex state is nonuniform in momentum space, but 



in real space transfers charge uniformly between layers. 
Our paper starts by describing a technical calculation 
that supports and elaborates on our prediction and then 
discusses anticipated properties of this state. 

Chiral two- dimensional electron system Hartree-Fock 
Theory — It is instructive to consider a class of chiral two- 
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FIG. 1: Pseudospin orientation in a graphene bilayer broken 
symmetry state. In this figure the arrows represent both the 
magnitude and the direction of the x-y projection of the pseu- 
dospin orientation n as obtained from a mean-field-theory cal- 
culation for a neutral, unbiased bilayer with coupling constant 
a — 1. The arrows are shorter in the core of the momentum 
space vortex because the pseudospins in the core have rotated 
spontaneously toward the z or — z direction. 



dimensional electron system (C2DES) models which in- 
cludes the continuum limits of single-layer and bilayer 
graphene sheets as special cases. These C2DES models 
have band Hamiltonians, 



H 



band — 



so(k c ) - 



cos(J^> fe ) t5 CT + sin(J0 fc ) T y a , 



V* 
2 



Cfc.o 



(1) 



2 



where a, a' are pseudospin labels, J is the chirality in- 
dex, r a is a Pauli matrix, fc c is the model's ultraviolet 
momentum cutoff, £o(fc c ) is the energy scale of the band 
Hamiltonian, and a sum over valley and spin components 
is implicit. In Eq. ([1]), V g is an external potential term 
which couples to the pseudospin magnet order parame- 
ter and corresponds in the case of bilayer graphene to 
an external potential difference between the layers. For 
single-layer graphene J = 1 and Eo(k c ) = hvpk c , while for 
bilayer graphene J — 2 and £o(fc c ) = fr 2 k 2 /(2m*). The 
dimensionless coupling constant of these C2DESs, which 
measures the interaction strength, can be defined as a = 
(e 2 k c /e)/£o(k c ) where e is the effective dielectric function 
due to screening external to the tt electron system. In the 
case of a single graphene layer a mono = e 2 /{ev-pK), while 
in the bilayer case, a^i = 2e 2 / (ev c h) , where v c — hk c /m*. 
If we choose Q hk c — ^/2m*^i for the bilayers, we have 
ot-H = a: m ono- Typically e ^ 2.5 which implies a dimen- 



sionless coupling constant a ~ 1. We use £o(fc c ) and k c 1 
as energy and length units in the rest of this paper. 

The C2DES Hartree-Fock Hamiltonian can be written 
(in dimensionless units) in the following physically trans- 
parent form: 



wm =- E 4,^'^v CT (fe)4,^, (2) 



j(0 



where B$ a = B^\k)5 a ,, a + B<*)(k) . Ta> 



A--K1 2tt \k-k'\ 2 



and the pseudospin field (fe) has band and interaction 
contributions, 



B®{k) = fc' 7 cos(J0 fc ) + a 



B«(fc) = fc J sin(J0 fc )+a 



Jly -\k-k'\d f(i) iy\ 



|fe'|< 



1 2tt |fc-fc'| 2 



2fc/ -|fc-fc'|d f(i) fUl 



d z k' e 



/diff \ K 



fe'|<i 2^ |fe-fc'| 2 



2 • J|fe'|<i 



d 2 fc' / 1 



2tt V fe - fc 



/i 



(4) 
(5) 
(6) 



Here i,j label the four valley and spin components of 
graphene's J = 1 and J = 2 C2DESs, rjW(fc) is the di- 
rection of B«(fc), /iuUfc') [/^(fe')] is the sum of (dif- 
ference between) low- and high-energy occupation num- 
bers, V g = V g /eo(k c ) is the gate potential in units of 
£o(fc c ), d = k c d is the distance between layers in units 
of k~ x in the bilayer case, and d = in the monolayer 
case. The term proportional to d on the right-hand side 
of Eq. ([6]) is the Hartree potential which opposes charge 
transfer between layers in the bilayer case. 

Local minima of the Hartree-Fock energy func- 
tional solve Eqs. ((4j- (J6j) self-consistently. Our fo- 
cus here is on the broken symmetry momentum- 
space vortex solutions in which pseudospins near 



k = tilt away from their band Hamiltonian x- 
y plane orientations toward the ±i direction, i.e., 
n«(fc) = (^ ) (fc)co8(J</> fe ),nW(fc)sm(J^ fe ),ni i) (A;)) 
with [n^(k)} 2 + [ni^ (fc)] 2 = 1. Pseudospin polariza- 
tion in the z-direction corresponds to charge transfer 
between layers. This ansatz yields effective magnetic 
fields whose x-y plane projections are parallel to the 
band Hamiltonian effective field. We find that fiW(fe) = 
{Bf (fc) cos( Jfa), Bf (fc) sin( J<j> k ), B[ l) (fc)) with 

Bf{k) = k J + a / dk'F ± (k,k')f^)n^^), (7) 
Jo 



where the exchange kernels in Eqs. (|7|) and j8]) are given by 



„ „ ,, / d<j) e~i d with q = q(k,k',<p) = Jk 2 + k' 2 -2kk' cos U ). The 

F ± (k,k) = k / cos(J0), 

Jo 2tt q 

F z (k,k') = k> ^ ^~ (9) 
Jo 2tt q 
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pseudospin-chirality induced frustration is represented by 
the factor cos(J^) in the first line of Eq. ([9]) which makes 
F± much smaller than F z . 

Pseudospin magnet phase diagram — We test the sta- 
bility of the "normal" state [n { z\k) = at V s = 0] 
solution of the Hartree-Fock equations by linearizing 
the self-consistency condition; nf 1 — B^rvf* /B± — > 
/ | n (i)_Q. This gives a fc-space integral equation 

n«(A) = V I dk' Mi j[k, k') n^(k'), (10) 

where 

MiA k,k')= ^ik^-kmf^) , 

fc J + a / dk"F ± (k,k") f^(k") 
Jo 

The normal state is stable when the largest eigenvalue 
of the linear integral operator M. in the right-hand side 
of Eq. (fT0|) is smaller than 1. Eigenvalues larger than 1 
are possible only because F± is smaller than F z , i.e., be- 
cause of pseudospin chirality. Phase diagrams for J = 2 
and J = 1 are plotted in Fig. [5] The pseudospin mag- 
net is more stable for larger coupling constant because 
it is driven by interactions, for larger J because the 
typical value of the band energy term proportional to 
k J decreases with J, and for smaller doping because 
/diff (^) i s then nonzero in a larger region of fc-space. The 
eigenvectors of M. specify the instability channel. The 
component-index structure of M implies that the eigen- 
values occur in groups of four, three of which [labeled an- 
tiferromagnetic (AF) in Fig. [2] correspond in bilayers to 
states with no net charge transfer, i.e., "Y^jfii (k) = 0. 
The ferromagnetic (F) instability in which all compo- 
nents are polarized in the same sense is opposed by the 
Hartree potential and delayed to larger coupling con- 
stant. For both AF and F instabilities, n z (k) is peaked 
at small k where the x — y pseudospin-planc exchange en- 
ergies are most strongly frustrated by chirality, and the 
kinetic energy term which opposes pseudospin magnetism 
is weakest. 

The physics that drives pseudospin magnetism in 
graphene bilayers is illustrated in Fig. [3] which partitions 
the condensation energy into band, Hartree, intralayer 
exchange, and interlayer exchange contributions. Spon- 
taneous layer polarization lowers the intralayer interac- 
tion energy at a cost in all other components. The over- 
all energy change is negative, and the broken symmetry 
state occurs, because the interlayer exchange energy of 
the normal state is weakened by the band-Hamiltonian 
induced frustration explained earlier. The cost in inter- 
layer exchange energy of pseudospin rotation is therefore 
much smaller than the gain in intralayer exchange en- 
ergy and the overall energy is reduced. The energy gain 
is considerably larger for AF broken symmetry states. 
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FIG. 2: (Color online) Phase diagram of C2DES's with J = 2 
and J — 1. For the J — 2 bilayer case we have taken d = 0.2. 
Pseudospin magnetism occurs at strong coupling a and weak 
doping /. (1 + / = n-f +n± where the pseudospin density n a = 
E M <4,i, CT Cfc,i, CT )/Af and N = E fe , s 1.) In the J = 2 bilayer 
case, the Hartree potential favors smaller total polarization so 
that the initial normal (N) state instability (blue separatrix) 
is to antiferromagnetic (AF) states in which the pseudospin 
polarizations of different valley and spin components cancel. 
At larger a, the normal state is unstable (green separatrix) to 
ferromagnetic (F) pseudospin states. In the J = 1 monolayer 
case d = so the phase boundaries (red separatrix) of F and 
AF broken-symmetry states coincide. 



Pseudospintronics — In Fig. [4] we illustrate typical re- 
sults for the pseudospin (layer) polarization £ = (nj — 
n l)/( n T + n i) °f a graphene bilayer as a function of gate 
voltage V g . The AF ground state at V g = 0, which has 
[Z 2 x SU{A))/[SU(2) x SU{2)] broken symmetry because 
of the freedom to choose any two spin or pseudospin com- 
ponents for (say) positive polarization, is gradually polar- 
ized by the gate voltage, but eventually becomes unstable 
in favor of polarizing more layers in the sense preferred 
by the gate voltage. At sufficiently strong gate voltages, 
the F ground state in which all layers are polarized in 
the same sense becomes the ground state. As the gate 
voltage is varied local minima of the Hartree-Fock energy 
functional become saddle points which are in the basin 
of attraction of another local minima. In this way, the 
self-consistent solutions exhibit hysteretic behavior. 

If only the fully polarized solutions existed these re- 
sults for pseudospin polarization as a function of gate 
voltage would be very much like the behavior expected 
for an easy-axis ferromagnet in an external magnetic 
field along the hard axis. In magnetic memories bista- 
bility enables information storage. In magnetic metal 
spintronics the dependence of the resistance of a circuit 
containing magnetic elements on the magnetization ori- 
entation of those elements gives rise to sudden changes 
in resistance with field (giant magnetoresistance) which 
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FIG. 3: (Color online) Condensation energy per electron 5e 
[in units of Eo(fcc)] as a function of a for an undoped (/ = 0) 
J = 2 C2DES. This figure shows results for both F (top) and 
AF (bottom) states. 



can be used to sense very small magnetic fields. Cur- 
rents running through such a circuit can also be used to 
change the magnetic state through spin-transfer torques. 
Pseudospin ferromagnetism in graphene bilayers could 
potentially lead to very appealing electrical analogs of 
both of these effects. Because of the collective behav- 
ior of many electrons, the pseudospin ferromagnet can 
be switched between metastable states with gate volt- 
ages that are much smaller than the thermal energy k^T, 
potentially enabling electronics which is very similar to 
standard complementary metal-oxide-semiconductor but 
uses much less power. This possibility is analogous to the 
property that a magnetic element can be switched be- 
tween magnetic states by Zeeman field changes that are 
extremely small compared to the thermal energy k-gT. 
Pseudospin-transfer torques, which are expected to oc- 
cur in electronic bilayer systems [l(| , can also be used to 
switch the pseudomagnetic state. 

Discussion — The proposals made here are based on ap- 
proximate calculations and must ultimately be confirmed 
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FIG. 4: (Color online) Metastable configurations of the pseu- 
dospin ferromagnet as a function of bias voltage V g [in units 
of eo(fcc)] with a = 1 and / = 0. We find self-consistent so- 
lutions of the gap equations {J}-® in which the pseudospin 
polarization has the same sense in all four components (fer- 
romagnetic), in three of the four components (ferrimagnetic), 
or in half of the four components (antiferromagnetic) . 



by experiment. Indeed, it is well known that Hartree- 
Fock theory (HFT) often overestimates the tendency to- 
ward broken symmetry states. For example HFT pre- 
dicts that a non-chiral 2DES is a (real-spin) ferromag- 
net at moderate coupling strengths, whereas experiments 
and accurate quantum Monte Carlo calculations suggest 
that ferromagnetism occurs only at a quite large value 
of the coupling constant [lj. Nilsson et al. have re- 
cently claimed that a similar ferromagnetic instability 
occurs in weakly-doped graphene bilayers, presumably 
only at a much stronger coupling constant than implied 
by HFT. We believe that the momentum-space vortex 
instability identified here, which is unique to the pecu- 
liar band-structure of bilayer graphene, is qualitatively 
more robust than the real-spin ferromagnetic instability. 
This should be especially true in neutral bilayers since 
the momentum-space vortex instability occurs at a cou- 
pling constant (a — > 0) for which correlation corrections 
to HFT are weak. This is not a strong-coupling insta- 
bility like ferromagnetism, but much more akin to the 
very robust attractive-interaction weak-coupling instabil- 
ity which leads to superconductivity. The condensation 
energy per electron associated with the formation of a 
momentum-space vortex core is ~ e 2 fc c /e, much larger 
than the ~ e 2 kp/e condensation energy for the spin- 
polarized state. Because this broken symmetry state is 
most robust for uniform neutral bilayers, the smooth but 
strong disorder potentials responsible for inhomogene- 
ity [12| in nearly neutral graphene sheets may need to 
be limited to allow this physics to emerge. 

Work at UT Austin was supported by the Welch Foun- 
dation, NRI-SWAN, ARO, and DOE. 
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